Remarkably Low KIR and HLA Diversity in Amerindians Reveals Signatures of Strong Purifying Selection Shaping the Centromeric KIR Region.

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival.

pMHC Structural Comparisons as a Pivotal Element to Detect and Validate T-Cell Targets for Vaccine Development and Immunotherapy-A New Methodological Proposal.

The search for epitopes that will effectively trigger an immune response remains the "El Dorado" for immunologists. The development of promising immunotherapeutic approaches requires the appropriate targets to elicit a proper immune response. Considering the high degree of HLA/TCR diversity, as well as the heterogeneity of viral and tumor proteins, this number will invariably be higher than ideal to test. It is known that the recognition of a peptide-MHC (pMHC) by the T-cell receptor is performed entirely in a structural fashion, where the atomic interactions of both structures, pMHC and TCR, dictate the fate of the process. However, epitopes with a similar composition of amino acids can produce dissimilar surfaces. Conversely, sequences with no conspicuous similarities can exhibit similar TCR interaction surfaces. In the last decade, our group developed a database and in silico structural methods to extract molecular fingerprints that trigger T-cell immune responses, mainly referring to physicochemical similarities, which could explain the immunogenic differences presented by different pMHC-I complexes. Here, we propose an immunoinformatic approach that considers a structural level of information, combined with an experimental technology that simulates the presentation of epitopes for a T cell, to improve vaccine production and immunotherapy efficacy.

Reconstructed Lost Native American Populations from Eastern Brazil Are Shaped by Differential Jê/Tupi Ancestry.

After the colonization of the Americas by Europeans and the consequent Trans-Atlantic Slave Trade, most Native American populations in eastern Brazil disappeared or went through an admixture process that configured a population composed of three main genetic components: the European, the sub-Saharan African, and the Native American. The study of the Native American genetic history is challenged by the lack of availability of genome-wide samples from Native American populations, the technical difficulties to develop ancient DNA studies, and the low proportions of the Native American component in the admixed Brazilian populations (on average 7%). We analyzed genome-wide data of 5,825 individuals from three locations of eastern Brazil: Salvador (North-East), Bambui (South-East), and Pelotas (South) and we reconstructed populations that emulate the Native American groups that were living in the 16th century around the sampling locations. This genetic reconstruction was performed after local ancestry analysis of the admixed Brazilian populations, through the rearrangement of the Native American haplotypes into reconstructed individuals with full Native American ancestry (51 reconstructed individuals in Salvador, 45 in Bambui, and 197 in Pelotas). We compared the reconstructed populations with nonadmixed Native American populations from other regions of Brazil through haplotype-based methods. Our results reveal a population structure shaped by the dichotomy of Tupi-/Jê-speaking ancestry related groups. We also show evidence of a decrease of the diversity of nonadmixed Native American groups after the European contact, in contrast with the reconstructed populations, suggesting a reservoir of the Native American genetic diversity within the admixed Brazilian population.

AVPR1b variation and the emergence of adaptive phenotypes in Platyrrhini primates.

Platyrrhini (New World monkeys, NWm) are a group of primates characterized by behavioral and reproductive traits that are otherwise uncommon among primates, including social monogamy, direct paternal care, and twin births. As a consequence, the study of Platyrrhine primates is an invaluable tool for the discovery of the genetic repertoire underlying these taxon-specific traits. Recently, high conservation of vasopressin (AVP) sequence, in contrast with high variability of oxytocin (OXT), has been described in NWm. AVP and OXT functions are possible due to interaction with their receptors: AVPR1a, AVPR1b, AVPR2, and OXTR; and the variability in this system is associated with the traits mentioned above. Understanding the variability in the receptors is thus fundamental to understand the function and evolution of the system as a whole. Here we describe the variability of AVPR1b coding region in 20 NWm species, which is well-known to influence behavioral traits such as aggression, anxiety, and stress control in placental mammals. Our results indicate that 4% of AVPR1b sites may be under positive selection and a significant number of sites under relaxed selective constraint. Considering the known role of AVPR1b, we suggest that some of the changes described here for the Platyrrhini may be a part of the genetic repertoire connected with the complex network of neuroendocrine mechanisms of AVP-OXT system in the modulation of the HPA axis. Thus, these changes may have promoted the emergence of social behaviors such as direct paternal care in socially monogamous species that are also characterized by small body size and twin births.

A new in silico approach to investigate molecular aspects of factor IX missense causative mutations and their impact on the hemophilia B severity.

Factor IX (encoded by F9) is a protein in the coagulation process, where its lack or deficiency leads to hemophilia B. This condition has been much less studied than hemophilia A, especially in Latin America. We analyzed the structural and functional impact of 54 missense mutations (18 reported by us previously, and 36 other mutations from the Factor IX database) through molecular modeling approaches. To accomplish this task, we examine the electrostatic patterns, hydrophobicity/hydrophilicity, disulfide, and H-bond differences of the Factor IX structures harboring the missense mutations found, correlating them with their clinical effects. The 54 mutated sequences were modeled and their physicochemical features were determined and used as input in clusterization tools. The electrostatic pattern seems to influence in disease severity, especially for mutations investigated in epidermal growth factors 1 and 2 (EGF1/2) domains. The combined use of all physicochemical information improved the clustering of structures associated to similar phenotypes, especially for mutations from GLA and EGF1-2 domains. The effect of mutations in the disease phenotype severity seems to be a complex interplay of molecular features, each one contributing to different impacts. This highlights that previous studies and tools analyzing individually single features for single mutations are missing elements that fulfill the whole picture.

A GWAS in Latin Americans highlights the convergent evolution of lighter skin pigmentation in Eurasia.

We report a genome-wide association scan in >6,000 Latin Americans for pigmentation of skin and eyes. We found eighteen signals of association at twelve genomic regions. These include one novel locus for skin pigmentation (in 10q26) and three novel loci for eye pigmentation (in 1q32, 20q13 and 22q12). We demonstrate the presence of multiple independent signals of association in the 11q14 and 15q13 regions (comprising the GRM5/TYR and HERC2/OCA2 genes, respectively) and several epistatic interactions among independently associated alleles. Strongest association with skin pigmentation at 19p13 was observed for an Y182H missense variant (common only in East Asians and Native Americans) in MFSD12, a gene recently associated with skin pigmentation in Africans. We show that the frequency of the derived allele at Y182H is significantly correlated with lower solar radiation intensity in East Asia and infer that MFSD12 was under selection in East Asians, probably after their split from Europeans.

Latin Americans show wide-spread Converso ancestry and imprint of local Native ancestry on physical appearance.

Historical records and genetic analyses indicate that Latin Americans trace their ancestry mainly to the intermixing (admixture) of Native Americans, Europeans and Sub-Saharan Africans. Using novel haplotype-based methods, here we infer sub-continental ancestry in over 6,500 Latin Americans and evaluate the impact of regional ancestry variation on physical appearance. We find that Native American ancestry components in Latin Americans correspond geographically to the present-day genetic structure of Native groups, and that sources of non-Native ancestry, and admixture timings, match documented migratory flows. We also detect South/East Mediterranean ancestry across Latin America, probably stemming mostly from the clandestine colonial migration of Christian converts of non-European origin (Conversos). Furthermore, we find that ancestry related to highland (Central Andean) versus lowland (Mapuche) Natives is associated with variation in facial features, particularly nose morphology, and detect significant differences in allele frequencies between these groups at loci previously associated with nose morphology in this sample.

Selection scan reveals three new loci related to high altitude adaptation in Native Andeans.

The Andean Altiplano has been occupied continuously since the late Pleistocene, ~12,000 years ago, which places the Andean natives as one of the most ancient populations living at high altitudes. In the present study, we analyzed genomic data from Native Americans living a long-time at Andean high altitude and at Amazonia and Mesoamerica lowland areas. We have identified three new candidate genes - SP100, DUOX2 and CLC - with evidence of positive selection for altitude adaptation in Andeans. These genes are involved in the TP53 pathway and are related to physiological routes important for high-altitude hypoxia response, such as those linked to increased angiogenesis, skeletal muscle adaptations, and immune functions at the fetus-maternal interface. Our results, combined with other studies, showed that Andeans have adapted to the Altiplano in different ways and using distinct molecular strategies as compared to those of other natives living at high altitudes.

Serotonin, behavior, and natural selection in New World monkeys.

Traits that undergo massive natural selection pressure, with multiple events of positive selection, are hard to find. Social behaviour, in social animals, is crucial for survival, and genetic networks involved in behaviour, such as those of serotonin (5-HT) and other neurotransmitters, must be the target of natural selection. Here, we used molecular analyses to search for signals of positive selection in the 5-HT system and found such signals in the M3-M4 intracellular domain of the 5-HT3A serotonin receptor subunit (HTR3A) in primates. We detected four amino acid sites with signs of putatively positive selection (398, 403, 432 and 416); the first three showed indications of being selected in New World monkeys (NWM, Platyrrhini), specifically in the Callitrichinae branch. Additionally, we searched for associations of these amino acid variants with social behavioural traits (i.e. sex-biased dispersal, dominance and social monogamy) using classical and Bayesian methods, and found statistically significant associations for unbiased sex dispersal (398L and 416S), unbiased sex dominance (416S) and social monogamy (416S), as well as significant positive correlation between female dispersal and 403G. Furthermore, we found putatively functional protein motifs determined by three selected sites, of which we highlight a ligand motif to GSK3 in the 416S variant, appearing only in Platyrrhini. 5-HT, 5-HT3A receptor and GSK3 are part of a network that participates in neurodevelopment and regulates behaviour, among other functions. We suggest that these genetic variations, together with those found in other neurotransmitter systems, must contribute to adaptive behaviours and consequently to fitness in NWMs.

Contrasting patterns of RUNX2 repeat variations are associated with palate shape in phyllostomid bats and New World primates.

Establishing the genetic basis that underlies craniofacial variability in natural populations is one of the main topics of evolutionary and developmental studies. One of the genes associated with mammal craniofacial variability is RUNX2, and in the present study we investigated the association between craniofacial length and width and RUNX2 across New World bats (Phyllostomidae) and primates (Catarrhini and Platyrrhini). Our results showed contrasting patterns of association between the glutamate/alanine ratios (Q/A ratio) and palate shape in these highly diverse groups. In phyllostomid bats, we found an association between shorter/broader faces and increase of the Q/A ratio. In New World monkeys (NWM) there was a positive correlation of increasing Q/A ratios to more elongated faces. Our findings reinforced the role of the Q/A ratio as a flexible genetic mechanism that would rapidly change the time of skull ossification throughout development. However, we propose a scenario in which the influence of this genetic adjustment system is indirect. The Q/A ratio would not lead to a specific phenotype, but throughout the history of a lineage, would act along with evolutionary constraints, as well as other genes, as a facilitator for adaptive morphological changes.

Molecular evolution and functional divergence of alcohol dehydrogenases in animals, fungi and plants.

Alcohol dehydrogenases belong to the large superfamily of medium-chain dehydrogenases/reductases, which occur throughout the biological world and are involved with many important metabolic routes. We considered the phylogeny of 190 ADH sequences of animals, fungi, and plants. Non-class III Caenorhabditis elegans ADHs were seen closely related to tetrameric fungal ADHs. ADH3 forms a sister group to amphibian, reptilian, avian and mammalian non-class III ADHs. In fishes, two main forms are identified: ADH1 and ADH3, whereas in amphibians there is a new ADH form (ADH8). ADH2 is found in Mammalia and Aves, and they formed a monophyletic group. Additionally, mammalian ADH4 seems to result from an ADH1 duplication, while in Fungi, ADH formed clusters based on types and genera. The plant ADH isoforms constitute a basal clade in relation to ADHs from animals. We identified amino acid residues responsible for functional divergence between ADH types in fungi, mammals, and fishes. In mammals, these differences occur mainly between ADH1/ADH4 and ADH3/ADH5, whereas functional divergence occurred in fungi between ADH1/ADH5, ADH5/ADH4, and ADH5/ADH3. In fishes, the forms also seem to be functionally divergent. The ADH family expansion exemplifies a neofunctionalization process where reiterative duplication events are related to new activities.

How strong was the bottleneck associated to the peopling of the Americas? New insights from multilocus sequence data.

In spite of many genetic studies that contributed for a deep knowledge about the peopling of the Americas, no consensus has emerged about important parameters such as the effective size of the Native Americans founder population. Previous estimates based on genomic datasets may have been biased by the use of admixed individuals from Latino populations, while other recent studies using samples from Native American individuals relied on approximated analytical approaches. In this study we use resequencing data for nine independent regions in a set of Native American and Siberian individuals and a full-likelihood approach based on isolation-with-migration scenarios accounting for recent flow between Asian and Native American populations. Our results suggest that, in agreement with previous studies, the effective size of the Native American population was small, most likely in the order of a few hundred individuals, with point estimates close to 250 individuals, even though credible intervals include a number as large as ~4,000 individuals. Recognizing the size of the genetic bottleneck during the peopling of the Americas is important for determining the extent of genetic markers needed to characterize Native American populations in genome-wide studies and to evaluate the adaptive potential of genetic variants in this population.

The Evolution of Science in a Latin-American Country: Genetics and Genomics in Brazil.

This article begins with a brief overview of the history of Brazil and that of Brazilian science, from the European discovery of the country in 1500 up to the early 21st century. The history of the fields of genetics and genomics, from the 1930s, is then first examined from the focal point of the lives and publications of the three persons who are generally considered to be the founders of genetics in Brazil (C. A. Krug, F. G. Brieger, and A. Dreyfus), and then by 12 other researchers up to 1999. The area of molecular genetics and genomics from 2000 to present is then described. Despite the problems of underdevelopment and the periodical political and economic crises that have affected life in Brazil, the fields of genetics and genomics in Brazil can be regarded as having developed at an appropriate pace, and have contributed in several major ways to world science.

Remembering the past - studies on evolution done by the genetics group at Universidade Federal do Rio Grande do Sul (UFRGS).

After a brief introduction about the factors that are involved in science development, and world and Brazilian evolutionary genetics, the studies developed in Porto Alegre in this area were reviewed. Four periods in the development of this group were distinguished: (a) Origins and first expansion (1949-1961); (b) Second expansion (1962-1988); (c) Third expansion (1989-2001); and (d) The last 15 years (2002-present). The international Porto Alegre Biological Evolution Workshops (PABEWs), with five biannual events from 2007 o 2015, were also mentioned. The final message stressed the importance of the maintenance of this and other Brazilian groups of research through adequate finance and recognition.

How strong was the bottleneck associated to the peopling of the Americas? New insights from multilocus sequence data

Abstract In spite of many genetic studies that contributed for a deep knowledge about the peopling of the Americas, no consensus has emerged about important parameters such as the effective size of the Native Americans founder population. Previous estimates based on genomic datasets may have been biased by the use of admixed individuals from Latino populations, while other recent studies using samples from Native American individuals relied on approximated analytical approaches. In this study we use resequencing data for nine independent regions in a set of Native American and Siberian individuals and a full-likelihood approach based on isolation-with-migration scenarios accounting for recent flow between Asian and Native American populations. Our results suggest that, in agreement with previous studies, the effective size of the Native American population was small, most likely in the order of a few hundred individuals, with point estimates close to 250 individuals, even though credible intervals include a number as large as ~4,000 individuals. Recognizing the size of the genetic bottleneck during the peopling of the Americas is important for determining the extent of genetic markers needed to characterize Native American populations in genome-wide studies and to evaluate the adaptive potential of genetic variants in this population.

Remembering the past - studies on evolution done by the genetics group at Universidade Federal do Rio Grande do Sul (UFRGS)

Abstract After a brief introduction about the factors that are involved in science development, and world and Brazilian evolutionary genetics, the studies developed in Porto Alegre in this area were reviewed. Four periods in the development of this group were distinguished: (a) Origins and first expansion (1949-1961); (b) Second expansion (1962-1988); (c) Third expansion (1989-2001); and (d) The last 15 years (2002-present). The international Porto Alegre Biological Evolution Workshops (PABEWs), with five biannual events from 2007 o 2015, were also mentioned. The final message stressed the importance of the maintenance of this and other Brazilian groups of research through adequate finance and recognition.

Molecular evolution and functional divergence of alcohol dehydrogenases in animals, fungi and plants

Abstract Alcohol dehydrogenases belong to the large superfamily of medium-chain dehydrogenases/reductases, which occur throughout the biological world and are involved with many important metabolic routes. We considered the phylogeny of 190 ADH sequences of animals, fungi, and plants. Non-class III Caenorhabditis elegans ADHs were seen closely related to tetrameric fungal ADHs. ADH3 forms a sister group to amphibian, reptilian, avian and mammalian non-class III ADHs. In fishes, two main forms are identified: ADH1 and ADH3, whereas in amphibians there is a new ADH form (ADH8). ADH2 is found in Mammalia and Aves, and they formed a monophyletic group. Additionally, mammalian ADH4 seems to result from an ADH1 duplication, while in Fungi, ADH formed clusters based on types and genera. The plant ADH isoforms constitute a basal clade in relation to ADHs from animals. We identified amino acid residues responsible for functional divergence between ADH types in fungi, mammals, and fishes. In mammals, these differences occur mainly between ADH1/ADH4 and ADH3/ADH5, whereas functional divergence occurred in fungi between ADH1/ADH5, ADH5/ADH4, and ADH5/ADH3. In fishes, the forms also seem to be functionally divergent. The ADH family expansion exemplifies a neofunctionalization process where reiterative duplication events are related to new activities.

Functional New World monkey oxytocin forms elicit an altered signaling profile and promotes parental care in rats.

The neurohormone oxytocin is a key player in the modulation of reproductive and social behavioral traits, such as parental care. Recently, a correlation between different forms of oxytocin and behavioral phenotypes has been described in the New World Monkeys (NWMs). Here, we demonstrate that, compared with the Leu8OXT found in most placental mammals, the Cebidae Pro8OXT and Saguinus Val3Pro8OXT taxon-specific variants act as equi-efficacious agonists for the Gq-dependent pathway but are weaker agonists for the ß-arrestin engagement and subsequent endocytosis toward the oxytocin receptor (OXTR). Upon interaction with the AVPR1a, Pro8OXT and the common Leu8OXT yielded similar signaling profiles, being equally efficacious on Gq and ß-arrestin, while Val3Pro8OXT showed reduced relative efficacy toward ß-arrestin. Intranasal treatment with either of the variants increased maternal behavior and also promoted unusual paternal care in rats, as measured by pup-retrieval tests. We therefore suggest that Val3Pro8OXT and Pro8OXT are functional variants, which might have been evolutionarily co-opted as an essential part of the adaptive genetic repertoire that allowed the emergence of taxon-specific complex social behaviors, such as intense parental care in the Cebidae and the genus Saguinus.

A tale of agriculturalists and hunter-gatherers: Exploring the thrifty genotype hypothesis in native South Americans.

OBJECTIVES: To determine genetic differences between agriculturalist and hunter-gatherer southern Native American populations for selected metabolism-related markers and to test whether Neel's thrifty genotype hypothesis (TGH) could explain the genetic patterns observed in these populations. MATERIALS AND METHODS: 375 Native South American individuals from 17 populations were genotyped using six markers (APOE rs429358 and rs7412; APOA2 rs5082; CD36 rs3211883; TCF7L2 rs11196205; and IGF2BP2 rs11705701). Additionally, APOE genotypes from 39 individuals were obtained from the literature. AMOVA, main effects, and gene-gene interaction tests were performed. RESULTS: We observed differences in allele distribution patterns between agriculturalists and hunter-gatherers for some markers. For instance, between-groups component of genetic variance (FCT ) for APOE rs429358 showed strong differences in allelic distributions between hunter-gatherers and agriculturalists (p = 0.00196). Gene-gene interaction analysis indicated that the APOE E4/CD36 TT and APOE E4/IGF2BP2 A carrier combinations occur at a higher frequency in hunter-gatherers, but this combination is not replicated in archaic (Neanderthal and Denisovan) and ancient (Anzick, Saqqaq, Ust-Ishim, Mal'ta) hunter-gatherer individuals. DISCUSSION: A complex scenario explains the observed frequencies of the tested markers in hunter-gatherers. Different factors, such as pleotropic alleles, rainforest selective pressures, and population dynamics, may be collectively shaping the observed genetic patterns. We conclude that although TGH seems a plausible hypothesis to explain part of the data, other factors may be important in our tested populations.